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SUMMARY 


Fall  out  contamination  deposited  on  the  roof  of  a  structure  is, 
in  many  cases,  the  source  of  the  primary  radiation  component  of  the  total 
dose  obtained  at  any  point  within  the  structure.  Experiments  have  been 
performed  in  which  the  doses  from  a  sou»ce  of  radiation  present  on  a  roof 
were  measured  in  many  locations  within  a  multi-story  building. 

This  report  presents  the  results  of  these  experiments  for  roof 
ana  floor  mass  thicknesses  of  48.6  and  97.  2  psf„  Comparisons  of  the 
experimentally  measured  gamma  doses  with  those  determined  theoretically 
have  been  shown  throughout  this  report.  Agreement  between  experiment 
and  theory  has,  in  general,  been  found  to  be  good. 


FOREWORD 


This  report  presents  results  of  an  experimental  evaluation  of  radia¬ 
tion  attenuation  and  distribution  caused  by  simulated  fallout  from  the  roof 
of  a  multi-story  structure  with  floors  of  varying  mass  thicknesses.  The  ex¬ 
periment  was  performed  during  the  period  August  1%5-1566  by  the 
CONESCO  Division  cf  Flow  Corporation  at  the  Protective  Structures  Dev¬ 
elopment  Center  (PSDC),  Fort  Belvoir,  Virginia. 

This  work  was  conducted  for  the  Office  of  Civil  Defense  through 
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is  Chief  of  its  PSDC  element.  All  Conesco  operations  under  this  contract 
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tal  Calibration,  and  Analysis  of  the  Radiation  Test  Facility  at  the  Protect¬ 
ive  Structures  Development  Center,  "  PSDC-TR-15,  "The  Barrier  Atten¬ 
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ral  Evaluation  of  Roof  Reduction  Factors.  " 
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NOMENCLATURE 


A  =  area  of  contaminated  field  (ft2) 

C  =  chafer-reader  cal  ibration  constant 

D  =  total  dose-roentgens 

D0  =  infinite  field  3  ft.  dose  rate 

e  =  eccentricity  of  structure  (W/L) 

H  =  height  of  source  above  ground  (ft.) 

I  =  specific  dose  rate  (R/hr)/(curie/?t2) 

L  =  building  length 

^(d/cosg)=  differential  angular  distribution  of  dose 

L(X)  =  the  total  detector  response  at  a  mass  thickness,  X,  from  an 
infinite,  plane  isotropic  source,  divided  by  the  total  de¬ 
tector  response  at  3  feet  in  air  from  the  same  source 

La  =  attenuation  introduced  by  a  barrier  adjacent  to  the  source 
plane 

=  attenuation  introduced  by  a  barrier  adjacent  to  the  detector 

Lc  =  attenuation  introduced  by  a  barrier  uniformly  distributed  be¬ 
tween  source  and  detector 

m  =  charger-reader  microamp  scale  reading 

n  =  2Z/L 

p  =  pressure  (in.  of  Hg) 

Rf  =  reduction  factor  =  (/ 

SQ  =  source  strength  in  curies 

S(X,cose)=  attenuation  of  parallel  monodirectional  radiation 
T  =  temperature  in  degrees  Rankine 

t  =  time  (hr) 

W  =  building  width 


effective  mass  thickness  of  barrier  (psf) 
source-to-detector  distance 

polar  angle  measured  relative  to  a  perpendicular  from  the 
detector  to  the  source  plane 

solid  angle  fraction 
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CHAPTER  1 


INTRODUCTION 


As  part  of  a  continuing  experimental  program  designed  to  verify 
the  theoretical  methods  utilized  in  computing  the  dose  rates  within  a 
structure  from  contamination  lying  on  the  roof,  the  CONESCO  Division 
of  the  Flow  Corporation  has  conducted  experiments  using  a  multi-story, 
concrete  structure.  The  results  of  these  experiments  and  their  comparison 
with  calculated  values  are  presented  in  this  report. 

The  initial  phase  of  the  roof  analysis  program  was  restricted  to 
experiments  performed  on  a  single-story  structure,  thus  concentrating  the 
barrier  near  the  source  of  radiation.  The  second  phase  of  the  program,  in 
which  floors  were  placed  within  the  structure,  provided  experimental  data 
which  were  compared  with  results  obtained  from  analytical  methods  ^  in 
which  the  attenuating  barrier  was  distributed  between  the  source  and  the 
detector. 


This  method  ,  in  its  current  form,  treats  barrier  masses  interposed 
between  roof  sources  and  detectors  as  being  uniformly  distributed  in  space. 
The  authors  of  this  method  have  assumed  that  this  treatment  is  more  repre¬ 
sentative  of  the  actual  situation  than  if  the  attenuating  mass  had  been  con¬ 
centrated  at  either  the  source  or  detector  locations. 


CHAPTER  2 

THEORETICAL  BACKGROUND 


2. 1  CALCULATIONAL  METHODS 

In  the  analysis  of  structural  shielding  from  nuclear  fallout, 
photon  attenuation  by  a  roof  is  presented  in  terms  of  a  reduction  factor. 

The  reduction  factor  is  defined  as  the  ratio  of  the  dose  at  any  location  from 
a  contaminated  area  of  arbitrary  geometry  to  the  dose  encountered  at  a 
point  3  feet  in  the  air  above  an  infinite,  smooth,  uniformly  contaminated 
plane  of  similar  density. 

Spencer^,  from  a  purely  theoretical  solution  of  the  transport 
equation  in  an  infinite  media,  has  developed  three  functions  representing 
different  assumptions  regarding  barrier  placement  that  are  to  be  employed 
as  approximations  to  roof  reduction  factors.  Each  of  these  functions  is  de¬ 
pendent  on  X,  the  barrier  mass  thickness,  and  o,  the  fraction  of  the  solid 
angle  subtended  at  the  detector  by  a  circular  source  field.  It  is  tacitly 
assumed  that  both  a  rectangular  plane  and  a  circular  source  plane  yield 
identical  reduction  factors  when  viewed  as  functions  of  a  solid  angle  frac¬ 
tion.  The  validity  of  this  assumption  in  the  case  of  the  barrier  concen¬ 
trated  near  the  source  has  been  adequately  demonstrated  experimentally 
for  the  case  of  a  roof  with  a  1.5  length-to-width  ratio  in  previous  work.3 

The  three  functional  expressions  developed  by  Spencer  are  pic¬ 
torial  ly  represented  in  Figure  2.  J.  Each  expression  represents  a  slightly 
different  geometric  situation.  A  detailed  analysis  of  the  functions  has 
previously  been  presented^'3  and  is  thus  only  outlined  briefly  here.  The 
expressions  are  utilized  to  represent  three  different  source-barrier-de¬ 
tector  configurations;  viz. ,  (I)  the  barrier  concentrated  near  the  source; 
(2)  the  barrier  concentrated  near  the  detector;  and  (3)  the  barrier  uniformiy 
distributed  between  source  and  detector.  The  functions  are  labelled,  res¬ 
pectively,  L(X)  L  (X,  u),  L(X)  (X,  u),  and  L(X)  Lc  (X,  u)  expressed  in 
terms  of  a  roof  reduction  factor.  The  term  L(X)  is  defined^  as  the  total  de¬ 
tector  response  beyond  a  barrier  thickness  X  from  an  infinite  plane,  iso¬ 
tropic  source  divided  by  the  total  detector  response  at  3  feet  in  the  air 
from  the  same  source. 


2.2  MASS  NEAR  SOURCE 


The  function  LQ(X,  u)  represents  the  reduction  in  detector  res¬ 
ponse  which  occurs  if  an  isotropic  detector,  separated  from  an  infinite 
plane  source  of  gamma  rays  by  a  barrier  thickness  X,  is  replaced  by  a  de¬ 
tector  responding  only  to  gamma  rays  incident  within  a  particular  cone  of 
directions  (Figure  2.  la).  Note  that  this  is  analogous  to  the  attenuating 
mass  being  concentrated  next  to  the  source,  because,  in  this  situation, 
the  attenuating  mass  defines  the  solid  angle.  The  function  is  expressed  in 
terms  of  the  differential  dose  angular  distribution  as  follows: 

rV 


L  a(X,u)  = 


eose^X'  COS0)  d  (COS0) 

$  (X,  cosq)  d  (c«e) 


where 


o  =  l  -  cos  q  is  the  solid  angle  fraction  of  the 
source  field  as  viewed  by  the  detector, 

j£  (X,  cos  q)  =  the  dose  angular  distribution  at  mass  thickness 

X  above  an  infinite  plane  source  of  contamination 
immersed  in  an  infinite  medium,  and 

d(cos  g)  =  the  differential  incremant  of  the  solid  angle. 


2.3  MASS  NEAR  DETECTOR 


A  second  estimation  of  the  dose  resulting  from  a  limited  area  of 
contamination  is  Lb(X, «).  This  function  expresses  a  fractional  reduction 
in  detector  response  which  occurs  if  an  infinite  plane  isotropic  source  is 
suddenly  constricted  to  emitting  radiation  only  into  a  limited  cone  of 
directions  about  the  perpendicular  toward  the  detector  (Figure  2.  lb).  This 
situation  is  approximately  equivalent  to  concentrating  the  barrier  at  the 
detector  as  photons  arriving  at  large  angles  find  it  increasingly  difficult  to 
penetrate  the  large  slant  distances  through  the  barrier.  This  function  is  ex¬ 
pressed  as: 
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Lb  (X,  u)  =  v 


JL  ^ 


cosq  )  d  (cos  e) 


SUK  COS 0  )  d  (cosq) 


where 


s(X,  cosq)  =  the  attenuation  provided  by  a  slab  of  thickness 
X  from  the  face  of  the  slab  which  is  irradiated 
by  parallel,  monodirectional  gamma  rays  of  in¬ 
cident  obliquity  q; 

v  =  a  proportionality  constant  so  that  integration 
over  all  obliquities  will  give  the  plane  source 
isotropic  source  results; 

j?<x,  cos  g)  -•  the  dose  angular  distribution  at  mass  thickness 
X  above  an  infinite  plane  source  of  contamina¬ 
tion  immersed  in  an  infinite  medium. 

2.4  ‘'ASS  DISTRIBUTED  BETWEEN  SOURCE  AND  DETECTOR 

The  expression  LC(X,  u)  is  the  ratio  of  the  response  of  an  isotropic 
detector  caused  by  a  circular  area  of  contamination  to  that  caused  by  an 
infinite  plane,  isotropic  source  of  gamma  radiation  (Figure  2. 1c).  Its 
analytical  expression  is  developed  from  the  equations  for  dose  above  an 
infinite  plane  field  of  contamination,  and  for  the  dose  above  the  same 
plane  containing  a  cleared  circular  area  centered  below  the  detector. 
These  equations  are  identical  and  a  function  only  of  the  slant  distance 
from  the  detector  to  the  edge  of  the  cleared  circle.  The  dose  rate  at 
height  X  (equivalent  mass  thickness)  above  an  infinite  field  of  contamin¬ 
ation  is  computed  as: 


■b 


cosq)  d  (cosq) 


with  quantities  as  previously  defined.  The  dose  rate  above  an  infinite 
plane  source  of  contamination  containing  a  cleared  circle  subtending 
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angle  q  below  the  detector  is  expressed  as: 

X  X 


L(- 


cosg' 


r\ 

L 

‘  (  x  \ 

1 

l^COS  Q']'  CC*  0 

d  (cos  q)  (2.4) 


-1 


By  subtraction/  the  dose  rate  from  a  circular  area  of  contamination  of  cone 
angle  o’  as  viewed  by  the  detector  is: 


Lc(X,o) 


m 


1 

i 

j 


(2.5) 


Roof  dose  measurements  in  a  multi-story  structure  have  Lsen  considered  to 
be  best  represented  by  the  L(X)  LC(X,  w)  function,  since  t.ies  barrier  is  gen¬ 
erally  located  in  finite  layers  between  the  source  and  detector.  This  func¬ 
tion  partially  accounts  for  the  radiation  from  an  isotropic  source  which 
has  been  scattered  to  the  detector  from  outside  the  solid  angle  subtended  by 
the  source  plane. 


2.5  CALCULATIONAL  ERRORS 

Spencer^  estimates  that  errors  in  the  calculations  caused  by 
truncation  of  the  momehts  solutions  and  the  correspondingly  required 
smoothing  of  the  angular  distribution  used  in  calculating  La(X/  u),  L^fX,  o), 
and  LC(X,  o)  range  between  0  and  25  percent.  On  the  other  hand,  m  the 
L(X)  calculations,  the  error  caused  by  moment  truncation  alone  reaches  a 
maximum  of  about  5  percent.  Errors  caused  by  inaccurate  input  informa¬ 
tion  (about  5  percent),  errors  inadvertently  introduced  when  drafting  the 
graphs,  and  those  caused  by  the  procedures  for  using  the  data  must  be 
added  to  these  estimates.  This  latter  error  is  traceable  to  the  inadequacies 
of  the  schematization  which  is,  in  part,  evaluated  in  this  study. 

2.6  SOLID  ANGLE  FRACTIONS 


As  previously  mentioned,  it  is  often  useful  to  represent  roof  re¬ 
duction  factors  as  functions  of  solid  angle  fractions.  The  solid  angle  frac¬ 
tion  is  defined  as  the  solid  angle  subtended  by  the  source  plane  at  the  de¬ 
tector  divided  by  2w.  For  a  circular  source  plane,  the  solid  angle  fraction 


6 


is  given  as: 


f 

j  d(cose)  =  1  -cose0  (2.6) 

where  0  is  the  angle  between  the  normal  to  the  center  of  the  disk  and  a 
line  drawn  from  the  detector  to  any  point  on  the  disk  (Figure  2.2). 

In  reality,  however,  most  roofs  are  best  represented  by  rectangular 
areas.  The  calculated  values  of  the  roof  reduction  factors  were 
developed  for  a  circular  geometry  and  the  assumption  is  implicitly  made 
that  the  detector  response  is  equivalent  from  a  circular  disk  source  and  a 
rectangular  source  subtending  the  same  solid  angle.  The  solid  angle  frac¬ 
tion  below  a  rectangle  is  represented  by 


2  1. 

u  = —  tan 

IT 


+  e^  +1 


(2.7) 


here 


e  =  ratio  of  length  to  width  of  rectangle 

n  =  perpendicular  distance  from  detector  to  rectangle  divided  by 
half  of  the  length  of  the  rectangle. 

2.7  PREVIOUo  WORK  IN  THIS  FIELD 

The  preceding  PSDC  report  on  roof  attenuation8  briefly  comment¬ 
ed  on  studies  previously  completed  on  this  subject.  These  include  work  by 
Eisenhauer,  Raso,  Clarke,6  Titus,  7  Batter,8  and  Schmoke9.  These 
experiments  were  conducted  in  single-story  structures  and  hence  are  not  di¬ 
rectly  comparable  with  this  work  except  on  the  third  floor  where  good  a- 
greement  was  achieved  with  the  L(X)  La  (X,y)  function.  The  experiment 
reported  herein  is  believed  to  be  the  first  detailed  study  conducted  on  a 
full  scale,  multi-story  building. 
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CHAPTER  3 


EXPERIMENTAL  ARRANGEMENT  AND  RESULTS 

3.1  BACKGROUND 

Experiments  to  determine  the  attenuation  of  gamma  rays  from  simu¬ 
lated  fallout  by  the  roof  of  a  structure  have  been  performed  at  the  Pro¬ 
tective  Structures  Development  Center.  Several  structural  geometries 
have  been  investigated  in  the  course  of  this  program  so  that  the  calcula- 
tional  methods  outlined  in  the  OCD  PM-lOO-l^  might  be  fully  evaluated. 

The  first  phase^  of  this  program,  was  primarily  concerned  with  single- 
story  structures.  In  that  phase  the  dose  rates  were  measured  below  a  source 
of  radiation  on  the  roof  of  a  box-type  structure  (consisting  of  four  walls  and 
a  roof  with  no  floors  or  partitions)  and  compared  with  existing  theory.  In 
the  second  phase  of  the  study  (described  in  this  report),  floors  were  intro¬ 
duced  in  the  structure,  thereby  producing  a  three-story  building.  Dose 
rates  were  then  measured  below  the  roof  source  on  each  of  the  stories.  The 
third  phase,  and  last  of  the  overall  roof  program,  will  retain  the  floors  in 
the  structure  and  introduce  interior  partitions  on  each  floor.  This  phase 
will  be  reported  separately. 

The  experimental  test  structure  and  the  method  of  obtaining  a  simu¬ 
lated  fallout  radiation  field  have  been  described  in  detail  in  a  previous  re¬ 
port^  and  hence  are  only  described  briefly  in  this  report. 

3.2  SOURCE  OF  RADIATION 

To  create  a  simulated  field  of  contamination  on  the  roof  of  the 
structure,  a  sealed  Co-60  source  was  pumped  at  a  constant  velocity  through 
a  length  of  polyethylene  tubing  pre-arranged  to  simulate  contamination  un¬ 
iformly  covering  the  particular  geometrical  configuration  of  interest.  In 
the  present  experiments,  two  geometries  were  investigated;  (1)  a  full  roof 
source,  i.e.,  a  rectangle  24  feet  wide  by  36  feet  long;  and  (2) a  circular 
disk  source  11,5  feet  in  diameter.  Figure  3. 1  illustrates  the  placement  of 
the  tubing  arrangement  on  the  roof  for  each  of  the  above  situations.  The 
smaller  area  disk  source  was  utilized  to  produce  test  results  for  smaller 
values  of  the  solid  angle  fraction  than  was  possible  with  the  full  rectangular 
roof  source. 


< 

u 

3.3  TEST  STRUCTURE 

V, 

The  experimental  test  structure  was  composed  of  a  steel  skeleton 
supporting  concrete  panels  as  walls  and  floors.  Figure  3.2  and  3.3  show 
the  test  structure  with  and  without  the  concrete  panels.  The  Interior  di¬ 
mensions  of  the  building  were  24  feet  In  width  by  36  feet  In  length  by 
36  feet  in  height.  The  bottom  face  of  the  floor  slabs  were  positioned  at 
11-,  23-,  and  35-fpot  elevations.  For  the  experiments  described  In  this  J 

report,  the  walls  of  the  test  structure  remained  at  a  constant  thickness  of 
4  inches.  Two  roof  and  floor  thicknesses  were  investigated;  viz.,  4  inches 
(48.6  psf)  and  8  inches  (97.2  psf). 

3.4  DETECTOR  POSITIONS 

In  the  experiments  in  which  the  disk  source  of  radiation  was  utilized 
on  the  roof,  measurements  of  the  dose  within  the  structure  were  made  only 
along  the  axis  intersecting  the  center  of  the  disk.  Since  the  disk  was  j 

positioned  in  the  center  of  the  rectangular  roof,  this  axis  corresponded  to  ! 

the  centerline  of  the  test  structure.  For  the  full  roof  source,  measurements 
were  made  not  only  along  this  centerline,  but  also  at  various  symmetric  off- 
center  locations.  A  floor  plan  of  the  test  structure  illustrating  the  center 
and  off-center  delimeter  locations  is  shown  In  Figure  3.4. 

3.5  ACCURACY  OF  EXPERIMENTAL  DATA 

To  determine  the  accuracy  of  the  data  obtained  from  the  roof  experi¬ 
ments,  the  errors  or  uncertainties  of  many  parameters  had  to  be  considered. 

Since  it  was  impractical  to  determine  experimentally  ail  of  the  uncertainties 
associated  with  weather,  exposure  time,  source  strength,  and  so  forth,  in  a 
completely  rigorous  way,  it  was  necessary  to  estimate  some  of  the  errors  and 
uncertainties  from  practical  experience,  A  detailed  analysis  of  those  errors 
is  presented  in  Appendix  B.  By  compounding  these  values  according  to  ac¬ 
cepted  principles,  the  standard  deviation  of  error  was  determined  and  is  pre¬ 
sented  in  Table  3.1. 
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Figure  3.2  Assembled  Test  Structure 
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Figure  3.3  Steel  Frame  of  Test  Structure 
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FIG.  3.4  -  SKETCH  OF  COORDINATE  CONVENTION  IN  TEST  STRUCTURE 


of  previously  published  data0  of  a  nearly  identical  situation,  with  and 
without  the  effect  of  the  supporting  beams,  indicated  however,  that  this 
effect,  while  large,  did  not  invalidate  the  comparison  of  experimental  re¬ 
sults  with  theoretical  ones  based  on  attenuation  through  concrete  alone. 


Tables  3.2  through  3.5  present  reduction  factors  measured  on  the 
roof  centerline  (perpendicular  to  the  plane  of  the  roof)  as  a  function  of 
the  detector  position  and  the  solid  angle  fraction  of  the  simulated  area  of 
contamination  (for  structures  with  both  4-  and  8-inch  floors)  for  source 
geometries  of  a  fully  contaminated  roof  (a  36  by  24  foot  rectangle)  and  an 
11.5-foot  diameter  circle.  The  calculated  reduction  factors,  L(X)Ld(X,  u) 
and  L(X)  Lc  (X,  u),  are  also  presented  in  these  fables,  together  with  the 
ratio  of  the  theoretical  to  experimental  reduction  factor  at  each  location. 
These  data  are  also  presented  graphically  in  Figures  3.5  and  3.6.  It 
should  be  noted  that  the  theoretical  values  of  reduction  factors  were  created 
by  summing  angular  distributions  that  were  computed  in  infinite  media  geo¬ 
metry.  Thus  both  functions  Lg  (X,  u)  and  LC(X,  u)  contained  radiation  back- 
scattered  from  the  atmospheric  mass  behind  the  roof  (skyshine).  The  ex¬ 
periments  described  in  this  study  did  not  contain  this  quantity  as  contam¬ 
ination  was  simulated  only  on  the  roof  of  the  structure,  and  the  majority 
of  the  skyshine  contribution  originates  a  mean  free  path  ('''-'450  feet  i  n 
air)  or  more  away.  An  estimate  of  this  quantity  however,  may  be  made  by 
the  following  equation  using  the  data  of  NBS  Monograph^  for  Co-60  rad¬ 
iation: 


Dsk  =  S(d)  Sa(d,u) 

where 

Dsk  =  skyshine  dose  from  an  infinite  field  of  Co-60 
radiation  normalized  such  that  the  total  dose 
at  a  height  of  3  feet  is  equal  to  1.0  R/hr 

d  =  mass  thickness  of  the  ceiling  expressed  in 
terms  of  equivalent  thickness  of  air=13.3 
times  thickness  in  psf 

u  =  solid  angle  fraction  through  which  the  de¬ 
tector  views  the  skyshine  (centered  about 
.  a  perpendicular  to  the  plane  of  contamination 

and  looking  away  from  it) 
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TABLE  3. 1 


EXPECTED  STANDARD  DEVIATION  OF  ERROR 
(Percent) 


Concrete  thickness 
Between 

Source  and  Detector 

Error  in  Dose  Rate 

Error  in  Reduction 
Factor 

inches 

PTF 

Percent 

Percent 

0 

0 

5.0 

6.2 

4 

48.6 

5.3 

6.5 

8 

97.2 

6.1 

7.1 

12 

145.8 

7.2 

8.1 

16 

194.4 

8.6 

9.4 

24 

243.0 

11.6 

12.2 

3.6  ANALYSIS  OF  EXPERIMENTAL  RESULTS 

In  this  section  the  results  of  the  experiments  described  above  are  pre¬ 
sented,  together  with  calculated  values  of  reduction  factors  obtained  by  a 
method  similar  to  the  OCD  PM-100-1  ^  using  attenuation  functions  repre¬ 
senting  the  mass  concentrated  at  the  source  and  evenly  distributed  between 
source  and  detector.  The  theoretical  values  of  these  functions  L(X)  Lq  (X,u) 
and  L(X)  LC(X,  to)  respectively,  were  taken  from  NBS  Monograph  42^  since 
it  was  desired  to  compare  the  experiments  with  a  Co-60  source  rather  than  a 
fallout  source  of  contamination.  This  was  done  so  that  conclusions  drawn 
were  not  affected  by  the  choice  of  energy  spectrum  and  hence  were  app¬ 
licable  to  fallout  as  well.  The  data  were  also  analyzed  and  compared  in 
the  light  of  previous  full-scale  experimental  work  performed  at  the  Pro¬ 
tective  Structures  Development  Center.  The  theoretical  function,  L(X)  Lb 
(X,  o)  is  not  investigated  in  this  report  since  the  physical  arrangement  used 
in  the  experiment  was  not  representative  of  this  source-barrier-detector 
configuration. 


3.6.1  Comparison  of  Experimental  Data  with  Theory 

An  "exact"  comparison  of  the  experimental  results  with  theory  was 
not  possible  because  of  the  effect  on  the  gamma  dose  rate  within  the  build¬ 
ing  of  the  non  uniform  mass  distribution  caused  by  the  steel  support  beams 
(not  "handleable"  in  the  theoretical  model)  in  the  test  structure.  Analysis 
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Reduction  Factor* 

Exp. 

L(X)La  IX,  u) 

L(X)  Lc(X,u) 

CONTAMINATION 


3.53x10"? 

5. 90x1 0"» 

4.50x10  , 

5.80x10  , 

4.61x10, 

5.65x10, 

4.57x10, 

5.50x10, 

4.59x10, 

5.40x10, 

4.33x10, 

5. 15x10, 

4.09x10  , 

4.90x10  , 

3.97*10, 

4.65x10, 

3.77x10  , 

4.45x10", 

3.60x10, 

4.20x10": 

3.45x10 

4. 00x1 O'* 

Standard  deviation  and  mean  value  of  ratlot 


9. 50x10": 
8.90x10": 


8.40x10  : 
7.90x10": 
7.50x10": 
7.ooxio": 

6.60x10“ 

6.30x10" 

6.00x10": 

5.70x10" 


Standard  deviation  and  mean  value  of  rafloj 


7. 84x1 0"? 

1.61x10"? 

9.53x10", 

1.52x10", 

1.00x10' 3 

1.44x10  , 

1.01x10~T 

1.39X10  , 

9.76x10"? 

9.37x10"? 

1.31x10  , 
1.28x10", 

8.90x10"? 

1.22x10, 

8.66x102? 

1.15x10  , 

8.41x10  ? 

1.11x10  , 

8. 79x1 O^? 

1.08x10  , 

7.74x10 

1.03xl0"j 

Standard  deviation  and  mean  value  of  mtlot 


2.05 

3.19 

1.59 

2.51 

1.44 

2.25 

1.38 

2.12 

1.34 

2.09 

1.37 

2.13 

1.37 

2.13 

1.33 

2.08 

1.32 

2.08 

1.23 

1.93 

1.33 

2.09 

1. 348+.  050 

2. 105+.060 
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TABLE  3,3 


Reduction  Factors  From  a  Circular  Roof  Source 
(48,6  psf  Roof  ond  Floors;  11.5'  dla.  source;  Centerline  Position) 


Detector  to 
Source  Plane 


Solid  Angle 
Fraction 


Reduction  Factors 


L(X)La  (X,  o)  L(X) Lc  (X,  u) 


Ratio  Theory  to  Exp, 


L(X)  La  (X 


Exp. 


ill 


*•**•••»**•******••*•  CONTAMINATION  ♦.»»»*»**»**•******* 


2.78x10  , 
3,59x10, 
3.41x10, 
3. 14x10, 
2.44x10, 
2.06x10, 
1.70x10, 
1.42x10": 


5.60x10  ; 
5.20x10"; 
4.60x10"; 
3. 90x10"; 
3.35x10" 
2.80x10"; 
2.40x10" 


2.00x10 


6.80x10 
6.30x10" 
5.40x10" 
4.50x10" 
3.80x10" 
3. 10x10* 
2.54x10" 
2. 10x10" 


Standard  deviation  and  mean  value  of  ratios  1.360+.060  1.470+;135 

////////////////  /  Four  Inches  Concrete  /  //  //  //  //"////  /  /  / 


0.075 

1.56x10"? 

2.90x10"? 

3.60x10^ 

1.86 

2.31 

0.066 

1.65x10", 

2.30x10  , 

3.10x10  , 

1.39 

1.88 

0.059 

1.52x10", 

2.03x10  : 

2.75x10  , 

1.34 

1.81 

0.053 

1.43x10", 

1.83x10  , 

2.49x10  , 

1.28 

1.74 

0.0475 

1.26x10, 

1.64x10  : 

2.20x10  , 

1.30 

1.75 

0.043 

1.19x10", 

1.50x10  : 

2.00x10  : 

1.26 

1.68 

0.039 

1.10x10, 

1.38x10  , 

1.81x10  : 

1.25 

1.65 

0.0354 

1.03x102! 

1.24x10  i 

1.67x10  : 

1.20 

1.62 

0.0324 

0.030 

9.66x10  \ 

9. 86x1 o"’ 

1.14x10  : 
1.07x10"° 

1.53x10": 

1.41x10"° 

1.18 

1.21 

1.58 

1.59 

Standard  deviation  and  mean  valuo  of  ratios 

1. 258+.  096 

1. 685+.  130 

LJJ-LLL. 


/  /  /  /  /  /  Four  Inches  Concrote  /  //  //  //  //  //  //  //  / 


1.54x10 

1.68x10" 

1.67x10" 

1.72x10" 

1.56x10" 

1.48x10" 

1.39x10" 

1.33x10" 

1.31x10" 

1.22x10" 

1.17x10" 


2.85x10 

2.68x10"' 

2.50x10"' 

2.31x10"' 

2.18x10“ 

2.04x10"' 

1.91x10“ 

1.81x10“ 

1.71x10“ 

1.62x10" 

1.55x10" 


4.55x10 

4. 22x10" 

3.94x10" 

3.62x10" 

3.40x10" 

3.20x10" 

3.00x10“ 

2.80x10" 

2.63x10" 

2.50x10" 

2.34x10" 


Standard  deviation  and  moon  valuo  of  ratios 


2. 125+.  145 
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TABU:  3.5 

Reduction  Factors  From  a  Circular  Roof  Source 
(97.2  psf  Roof  and  Floors;  1T.51  dia.  source;  Centerline  Position) 
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FIG.  3,5  -  Calculated  and  Experimental  Reduction  Factor*  for  48.6  p»f  and  97,2  p«f 
Mom  Thickneu  Floor*.  Full  Roof  Contamination 
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- L(X)L0<X,u) 

X  X  XXX  XX Experimental  97.2  p»f 
AAAAA  Experimental  48.2  p*f 


FIG.  3.6  -  Calculated  and  Experimental  Reduction  Factor*  for  48.6  p*f  and 
pif  Floor*,  11.5  ft  Diameter  Circular  Area  of  Contamination 


Sa(d, u)  =  angular  distribution  of  skyshine  dose 
(note  Sa  (d,  u)  =  1  at  u=l) 

If  Djk  Is  evaluated,  it  is  found  to  be  very  small  in  all  cases  covered 
in  this  study  (thickness  >  48.6  psf).  The  maximum  effect  expressed  in 
terms  of  equivalent  massthickness  of  the  roof  was  of  the  order  1  psf  for  the 
case  of  the  48.6  psf  roof.  Thus,  the  lack  of  adequate  simulation  of  sky  - 
shine  in  the  experiment  caused  the  roof  to  appear  to  be  a  maximum  of  a  - 
bout  1  psf  more  in  thickness  for  the  thin  roof  (48.6  psf  total  thickness)  and 
the  apparent  increase  in  thickness  was  negligible  in  all  other  cases. 

Tables  3.2  through  3.5  indicate  thrrt  either  of  the  theoretical  functions 
always  yields  a  conservative  value  (calculated  reduction  factor  greater  nu¬ 
merically  than  measured)  of  the  roof  reduction  factor.  The  attenuation  af¬ 
forded  by  the  support  beams  of  the  structure  could,  however,  account  fot 
at  least  part  of  this  difference.  It  Is  important  to  determine  if  the  theore¬ 
tical  value  of  L(X)  Lc  (X,  u)  or  L(X)  La  (X,  o)  would  or  would  not  under¬ 
estimate  the  roof  reduction  factors  within  the  structure  if  the  support  beams 
had  been  absent. 

The  effectveness  of  any  given  amount  of  mass  in  attenuating  radiation 
emanating  from  an  isotropic  plane  source  is  decreased  if  the  mass  is  lumpfed. 
in  discrete  masses  rather  than  spread  uniformly  to  attenuate  all  radiation. 

This  is  because  radiation  is  attenuated  exponentially,  while  mass  is  distri¬ 
buted  arithmetically.  For  example,  the  mass  required  to  reduce  the  dose 
rate  eminating  from  a  uniformly  contaminated  infinite  place  of  1.25  Mev 
radiation  by  a  factor  of  ten  is  41  PSF.  If  this  same  total  mass  per  unit  area 
were  lumped  into  areas  of  zero  PSF  and  410  PSF  thickness,  the  percentage 
of  the  contaminated  plane  that  was  shielded  by  410  PSF  would  by  only  10%, 
thus  the  total  dose  rate  measured  above  the  plane  with  lumped  masses  would 
by  slightly  more  than  90%  of  that  from  the  unshielded  source  plane  or  nine 
times  greater  than  if  the  same  amount  of  mass  had  been  uniformly  distributed 
over  the  source  plane. 

Thus  a  conservative  estimate  of  the  worth  of  the  steel  beams  could  be 
gained  by  "smearing"  the  steel  over  the  entire  area  of  the  floor.  In  other 
words,  the  total  amount  of  steel  per  floor  is  calculated  and  the  smeared 
thickness  of  an  equal  weight  steel  plate  is  obtained.  This  must  then  by 
multiplied  by  0, 931  to  correct  for  the  deficiency  of  electrons  in  iron  per 
unit  mass  compared  with  concrete^.  This  smeared  thickness  has  been  cal¬ 
culated  to  be  approximately  9.6  psf  per  floor. 
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Table  3.6  summarizes  the  data  of  Tables  3.2  through  3.5  in  the  form 
of  mean  values  and  standard  deviations  of  the  ratio  of  a  theoretically  cal¬ 
culated  dose  to  that  measured  experimentally.  An  effective  floor  thick¬ 
ness  may  thus  be  determined  for  each  method  (that  required  to  reduce  the 
mean  value  of  the  ratio  of  calculated-to-measured-dose  rate  to  1.0)  of 
computation.  Comparable  cases  are  presented  in  Table  3.6  together  with 
the  difference  between  this  "effective"  floor  thickness  and  the  actual  thick¬ 
ness  of  concrete.  Inspection  of  these  differences  indicates  that  even  if  the 
support  beams  were  as  effective  as  if  their  mass  were  smeared  over  the  entire 
floor  areas/  the  function  L(X)  LC(X;  u)  would  serve  as  a  maximum  upper  (imit 
to  the  roof  reduction  factors.  However,  since  the  effect  of  inhomogene  ous 
mass  distribution  must  be  less  than  if  the  steel  were  smeared,  L(X)La  (X,  u), 
in  all  probability,  represents  a  realistic  value. 

O 

Previous  work0  in  idealized  geometry  with  the  mass  concentrated  ad¬ 
jacent  to  the  source  (and  with  a  scattering  mass  behind  the  source  to  repro¬ 
duce  the  skyshine  effect)  indicated  that  the  best  agreement  between  theory 
and  experiment  could  be  obtained  in  the  case  of  a  48.6  psf  (4-inch  concrete) 
shield  if  theoretical  calculations  were  based  on  55  psf;  similarly  for  a  97.2 
psf  shield  calculation  for  100  psf  and  for  a  145.8  chield  calculation  for 
146  psf.  If  these  higher  values  are  taken  as  the  effective  thickness  of  the 
concrete  shield,  an  estimate  of  the  steel  effectineness  may  be  made.  This 
effect  of  the  steel  is  presented  in  the  right  hand  column  of  Table  3.6.  Note 
that  in  ihe  previous  work  the  mass  was  concentrated  only  at  the  source,  thus 
this  is  tfwonly  situation  that  may  be  compared  directly.  Since  the  difference 
between  experimentally  measured  dose  rates  and  those  theoreticaily  calcul¬ 
ated  (expressed  as  equivalent  mass  thickness  in  Table  3.6)  is  greater  than  the 
effect  of  the  support  beams,  with  the  experimentally  measured  dose  always 
lower  titan  that  calculated,  both  of  the  methods  of  calculation  (L(X)La(X,  u) 
and  L(X)  Lc  (X,  u)  yield  conservative  results. 

The  function  L(X)  Lq  (X,  u),  however,  seems  more  representative  of 
the  actual  reduction  factor  than  that  presently  used^  L(X)  Lc  (X,  u).  Since 
the  use  of  L(X)  La  (X,u)  may  increase  predicted  protection  factors  signif¬ 
icantly  in  some  cases,  further  experimental  and  theoretical  work  on  inho- 
mogenous  roof  slabs  is  warranted. 

The  data  of  this  experiment  may  also  be  used  to  illustrate  the  differ¬ 
ence  represented  by  the  theoretical  functions  L(X)  Lc  (X,  o)  and  L(X)L<j(X,  w). 
These  functions  L(X)  Lc  (X,  u)  and  L(X)  La  (X,  o)  represent,  respectively, 
the  case  in  which  the  barrier  is  uniformly  distributed  between  source  and  de- 
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tector  and  the  case  in  which  the  barrier  is  concentrated  near  the  source. 

The  former  case  predicts  a  higher  roof  contribution  than  the  latter  over 
the  range  of  the  investigation  where  X  varied  from  48.6  psf  to  291.6  psf. 
Experimental  verification  of  this  effect  is  presented  for  8  inches  of  concrete 
in  Figure  3. 7  where  measured  values  of  the  roof  reduction  factor  for  the 
case  of  two  4-inch  floors  separated  by  12  feet  of  open  space  and  the  case 
of  an  8-inch  solid  roof  adjacent  to  the  source  in  an  open  structure  of  simple 
box  geometry  are  presented.  The  difference  between  these  two  sets  of  ex¬ 
perimental  data  is  observed  to  be  similar  but  not  so  great  as  theory  would 
predict.  Thus,  here  again  there  exists  an  indication  that  the  L(X)  Lc  (X,u) 
function  may  be  overly  conservative  in  predicting  the  effect  of  discrete 
floor  masses. 

3.6,2  Off-center  positions 

In  the  case  of  the  entire  roof  of  the  structure  being  covered  with 
simulated  contamination,  detectors  were  placed  along  both  the  vertical 
centerline  and  at  various  off-center  locations.  Figure  3.4  shows  the  off- 
center  dosimeter  locations  in  the  test  structure.  The  data  from  these  ex¬ 
periments  expressed  as  roof  reduction  factors,  for  the  off-center  locations 
are  presented  in  Table  3.7  for  both  structures  (4-  and  8-inch  floor  thick¬ 
nesses)  tested.  A  typical  comparison  of  experiment  with  theory  is  shown 
in  Figure  3.8.  Where  the  theoretical  values  calculated  by  the  "fic¬ 
titious  building"  method  described  in  Append  No  new  information 
is  evidenced  in  this  comparison;  however,  it  4ing  to  note  that  ex¬ 

periment  and  theory  compare  in  the  same  fashiu  >rf-center  positions  as 
for  those  on  the  vertical  centerline.  This  observation  was  previously  made 
in  Ref.  3  in  which  only  the  case  of  the  mass  adjacent  to  the  source  (i.e., 
single-story  structures)  was  investigated.  The  conclusion  previously  made 
—that  the  "fictitious  building"  method  of  calculating  off-center  roof 
reduction  factors  is  at  least  as  valid  as  the  method  used  for  centerline  po¬ 
sitions - remains  valid. 

Diagonal  traverses  of  the  roof  reduction  factors  for  the  full  roof  source 
for  various  heights  in  each  structure  are  also  presented  in  Figure  3. 8.  The 
experimental  points  are  shown  together,  with  the  calculated  values.  Here 
again  the  same  conclusion  as  above  can  be  drawn  from  these  plots  as  to 
the  accuracy  of  the  "fictitious  building"  method. 
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-FIG.  3. 7  -  Comparison  of  Theory  with  Experiment  for  97.2  psf  Floor  Thickness 
Between  Source  and  Detector  for  Both  a  Single  and  Divided  Slab 
(Two  48.6  psf  Slabs  Separated  by  a  Distance  of  12  ft) 
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Typical  Roof  Reduction  Factors  Off  Center  Positions 
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FIG.  3,8  -  Comparison  of  Calculated  and  Experimentally  Measured  Roof  Reduction  Factors 

Off-Center  Positions 


CHAPTER  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  study  was  to  determine  what,  if  any,  significant 
differences  exist  in  the  dose  rate  from  roof-based  sources  of  contamination 
when  the  attenuating  mass  was  distributed  in  finite  layers  between  source 
and  detector  rather  than  concentrated  at  the  source  and  to  compare  the 
dose  rates  determined  with  those  predicted  theoretically.  The  interpre¬ 
tation  of  the  results  obtained  is,  unfortunately,  somewhat  clouded  by  the 
seeming  I  y  minor  inhomogeneity  of  the  ceiling  slabs  and  steel  supports. 

This  disparity,  caused  by  the  support  structure  was  of  the  same  order  of 
magnitude  as  the  effect  that  was  measured.  Conclusions  as  to  the  validity 
of  currently  used  theoretical  estimates,  however,  were  made  by  interpret- 
i  ng  the  results  obtained  in  this  study  in  the  light  of  previously  published 
work.  ^ 

4.1  CONCLUSIONS 

The  following  conclusions  resulted  from  this  study. 

1.  The  experimentally  measured  roof  reduction  factors  illustrate 
the  same  dependence  on  the  solid  angle  of  the  contaminated  areas  of 
those  computed  with  either  function  L(X)LC  (X,  u)  or  L(X),  La  (X,  u). 
Measured  values  were  always  significantly  lower  than  the  ones  pre¬ 
dicted.  When  the  effect  of  the  floor  support  beams  was  neglected, 
measured  values  were  as  much  as  25  and  50  percent  lower,  respectively, 
than  those  computed  using  L(X)  LQ(X,u)  and  L(X)  LC(X,  u). 

2.  The  current  use  of  the  function  L(X)  LC(X,  o)  to  represent 
the  roof  reduction  factor  seems  overly  conservative.  The  function 
L(X)  Lq  (X,  u)  presents  a  conservative  estimate  of  the  dose  rates 
measured. 

3.  Off-center  roof  reduction  factors  calculated  by  the  "fictitious 
building"  method  portray  the  same  relative  agreement  with  experiment 
as  do  the  centerline.  This  indicates  that  the  off-center  reduction  fac¬ 
tor  calculations  are  as  valid  as  these  on  the  structure  centerline. 
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4.2  RECOMMENDATIONS 

The  following  recommendations  are  offered  as  a  result  of  this 
study. 

1.  A  series  of  experiments  should  be  undertaken  upon  roofs  of 
controlled  non-uniform  thickness  to  ascertain  the  best  method  of  cal¬ 
culating  the  effect^  the  inhomogeneous  distribution  of  mass  in  a 
structure  roof  since  it  has  been  shown  to  produce  fairly  major  discre¬ 
pancies  in  calculated  reduction  factors. 

2.  Further  measurements  should  be  made  of  the  attenuation 
afforded  by  roof  and  ceiling  slabs  for  idealized  geometry  in  tho  ab¬ 
sence  of  mass  discontinuities, 

3.  Unless  further  experiments  indicate  differently,  it  is  tenta¬ 
tively  recommended  that  the  function  L(X)  La(X,u)  rather  than 
L(X)  L  (X,o)  be  used  to  represent  a  more  realistic  estimate  of  the 
dose  attributed  to  contamination  existing  on  the  roof  of  a  structure. 
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APPENDIX  A 
OFF-CENTER  POSITIONS 


Off-center  reduction  factors  resulting  from  roof  sources  can  be  cal¬ 
culated  using  the  "fictitious  building"  method.  This  method  consists  of 
dividing  the  roof  of  the  structure  being  considered  into  two  or  more 
imaginary  roofs  that  are  symmetrical  about  the  point  in  question,  and 
hence  calculable.  It  is  then  assumed  that  the  average  of  the  imaginary 
roof  reduction  factors  results  in  the  reduction  factor  of  the  off-center  point 
considered. 

As  an  example  of  this  procedure,  consider  a  point  located  at  position 
(6, 9)  (x,  y  coordinates  from  center  of  structure)  below  a  rectangular  24-  by 
36-foot  roof  of  145. 8  psf  thickness.  The  roof  area  must  be  divided  into  four 
rectangles,  each  of  which  is  symmetrical  about  the  point  of  interest.  Each 
of  these  rectangles,  labelled  Fictitious  Buildings  A,  B,  C,  and  D  in  Fig. 

A-l,  represents  four  times  the  dose  that  would  be  received  by  a  detector 
from  the  upper  left,  upper  right,  lower  left,  and  lower  right  quadrants, 
respectively,  if  the  coordinate  axis  were  centered  upon  the  point  of  interest. 
The  calculation  of  the  dose  from  each  of  these  centered  rectangles  then 
proceeds  in  the  manner  prescribed  by  the  OCD  PM  100-1  and  is  presented 
below: 
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The  total  of  the  reduction  factor  column  (L(X)La(X,  u)  divided  by  4  thus 
represents  the  expected  reduction  factor  below  the  point  6,  9.  Thus: 

Total  L(X)La(X,o)  =  0.0121 


Reduction  factor  =  -  =  0.U030 


J  Design  and  Review  of  Structures  for  Protection  from  Fallout  Gamma  Rad 
iation,  OCD  PM-100-1  (February  1 965) 
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APPENDIX  B 


ERROR  ANALYSIS 


A.  GENERAL 

The  accuracy  of  the  experimental  data  presented  in  this  report  is 
determined  by  the  errors  (or  uncertainties)  associated  with  each  of  the  para¬ 
meters  involved  in  the  determination  of  an  experimental  value.  These 
errors  are  determined  individually  and  then  compounded  according  to  the 
accepted  statistical  principles  to  ascertain  an  estimate  of  the  total  error  to 
be  associated  with  the  experiment. 

B.  STANDARD  DEVIATION  IN  SPECIFIC  DOSE  RATE 


The  specific  dose  rate  at  any  point  is  determined  by  the  following 

equation: 

D  =  -ML  (B-l) 

10J  pest 

where 

A 

D  =  specific  dose  rate  (R/hr/curie/fT) 

M  "  scale  reading  on  the  dosimeter  charger-reader  (pa) 

A  =  area  of  simulated  source  field  (ft^) 

T  =  temperature  (°R) 

p  =  barometric  pressure  (in.  Hg) 

c  =  charger-reader  cal  ibration  constant 

t  =  exposure  time  (hours) 

s  =  source  strength  (curies) 

The  error  in  any  measurement  of  D  is  related  to  the  errors  in  the 
components  by  the  total  differential  of  D  or 


£D  = 


d  P.  AM  +  -4-?  M  +-L.P  A?  +— 

TTvT  ™  TA  8P^  8 


95 


or 


(P-2) 
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By  evaluating  the  partial  derivatives,  Eq.  (B-2)  is  equivalent  to: 


AD  =  D 
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M  AT 

"7T+  ~T 


-  £p_ 
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w 

(B-3) 


The  expected  standard  deviation  (cq)  of  an)  ?  many  measurements  of 
D  is  related  to  the  standard  deviations  of  its  ent  parts  by 


When  all  of  the  panmeters  on  the  right  side  are  independent  of 
one  another,  it  is  often  more  convenient  to  speak  of  the  error  in  terms  of  a 
fractional  standard  deviation. 

Thus: 

ffD  =  D 

^M 

(B-5) 


where  is  the  mean  reading  of  D 

Generally  D/Dm  should  be  small  and  may  be  assumed  approximately  equal 
to  1. 

C.  MEASUREMENT  ERRORS. 

The  errors  in  the  individual  variables  which  together  yield  a 
specific  dose  rate  have  been  measured  or  estimated  as  follows:  (Reference 
1) 


1.  An  Experimental  Evaluation  of  Roof  Reduction  r'actors,  by  C.  H. 
McDonnell,  J.  Velletri,  PSDC-TR-16 
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1.  Microamp  Reading  -  The  standard  deviation  of  the  microamp  reading  on 
the  charger-reader  has  been  found  to  consist  of  two  components: 

a)  Scale  reading  error:  0. 5  pA 

b)  Normal  measured  dosimeter  group  variation  (see  Appendix  C): 

1.0  pa. 

Thus  Cw  is  estimated  to  be  approximately  1.0 pa  as  the  scale  reading 
error  is  included  in  this  dosimeter  variation. 

2.  Simulated  Source  Area  -  Extreme  care  was  taken  in  laying  out  the 
tubing  to  simulate  radiation  fields.  The  error  is  thus  believed  to  be 
quite  small.  A  value  of  of  0.001  (0. 1  percent)  has  been  assumed 
as  typical. 

3.  Absolute  Temperature  -  The  largest  source  of  error  in  determing  a 
temperature  at  whicli  an  experiment  was  run  is  due  to  the  variation  of 
temperatuje  with  time  during  the  course  of  an  exposure.  The  total 
error  in  temperature,  consists  of  two  components,  scale  reading  error 
and  time  variation  of  temperature  during  the  experiment.  The  stand¬ 
ard  deviation  of  this  value,  a-j-  is  estimated  to  be  approximately  3°R. 

4.  Barometric  Pressure  -  Similarly,  a  small  error  is  introduced  by  the  time 
variation  and  scale  reading  error  of  the  barometric  pressure.  These  are 
assumed  to  be  of  approximately  equal  magnitude.  The  standard  de¬ 
viation  a  is  thus  estimated  to  be  0.014  in.  Hg. 

r 

5.  Charger-Reading  Calibration  Constant  -  The  charger-readers  used  in 
the  experiments  are  calibrated  as  carefully  and  often  as  possible. 
Certain  errttrs,  however,  are  inherent  in  the  calibration  of  the  in¬ 
struments  and  the  calculation  of  the  calibration  constant.  Eh'  'lose 
observation  of  the  results  of  the  calibration  of  various  charger-readers 
with  various  types  of  dosimeters,  it  has  been  found  that  ac/c  ‘s  a.P" 
froximately  equal  to  0. 035  or  3. 5  percent,, 

6.  Source  Strength  -  All  the  sources  were  calibrated  using  NBS-calibrated 
Victoreen  chambers  and  charger-readers.  In  addition  to  small  errors  in 
these  instruments,  there  are  also  errors  introduced  by  the  measurement 
of  source-to-dfetector  distances  and  air  density  determinations.  It  has 
been  estimated  that  an  error  of  2  percent  (  oj  s  =  0. 02)  is  to  be  as  - 
soa'atedwith  the  source  strength  measurement. 

7k  Exposure  Time  -  Measurement  of  the  time  in  which  the  source  is  ex- 
posed  is  done  with  a  stopwatch.  The  error  induced  is  primarily  that  of 
reaction.  For  any  given  run,  ct  is  estimated  to  be  1.0  second  or 
0. 00028  hour. 


36 


D,  SUMMARY  OF  INDIVIDUAL  ERRORS 


A  summary  of  the  individual  errors  follows: 


Microamp  reading: 

aM  = 

1  ua 

Source  area: 

Q 

> 

> 

II 

0.001 

Absolute  temperature: 

/ 

°1  = 

3°R 

Barometric  pressure: 

a  = 

p 

0.014  in.  Hg. 

Calibration  constant: 

a/c  = 

0.035 

Source  strength: 

’/•  = 

0.02 

Exposure  time: 

0. 00028  hour. 

Thus,  for  any  experiment/  equation  (B-5)  reduces  to: 


These  errors  are  the  rame  for  all  experimental  values  throughout  this  re¬ 
port.  It  is  felt  that  they  accurately  reflect  the  range  within  which  at 
least  68. 2  percent  of  the  average  value  lie. 


E.  ERROR  ASSOCIATED  WITH  SPECIFIC  DOSE  RATE  FROM  UNATTEN¬ 
UATED  ROOF  SOURCES 

With  each  experimental  exposure  to  determine  dose  rates  within  the 
structure,  measurements  were  also  made  above  the  structure.  This  proced¬ 
ure  provided  unattenuated  dose  rates  at  various  source-to-detector  dis¬ 
tances  on  the  axis  defining  the  centerline  of  the  source  field.  Many  vari¬ 
ations  were  necessarily  noticed  in  the  individual  parameters  of  Eq.  (B-6). 
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The  error  measurement  described  above  was  performed  on  many  of  these 
points  and  the  following  are  offered  as  examples  of  the  use  of  the  method 
of  calculating  aQ 


1.  Full  Roof  Source 


D 


M 


2  = 

8  feet  above 

source 

M  = 

50  pA 

1  =  0.02 
rr 

T  = 

510  °R 

1  Yi 

=  0. 00588 

P  = 

29. 78  in.  Hg 

°*014  -  0.00047 

P 

t  = 

0. 27  hour 

0.00028  =  0>001037 

D  = 

77.  93  R/fhr/ < 

curie/  ft^j 

dm  = 

78.5  R^hr/ curie/ ft2] 

770  93 
75 75~ 

r  2  2  2 

[0. 001626  +  (.  02)  +  (0. 00588)  +  (0. 00047) 

(0.001 03 7)2 


1/2 


M 


=  0. 0451  =  4. 51  per  cent 


2.  Diameter  Circle  (11.5  feet) 

Z  =  8  feet  above  source 

M  =  38  p  a 

T  =  508°R 

p  =  29. 55  in.  Hg 


1 

77T 

3_ 

T 

0.014  = 
P 


0. 0263 
0. 0059 
0.00045 
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*? 


t  =  0.0897  hr  0.00028  =  0.00312 

t 

D  =  21.63  hr/  curie/  ft^ 

=  21,6  hr/  curie/ 


_  21.63 
"7T5" 


=  0. 0487  =  4. 87  percent 


The  mean  dose  rate,  D^,  was  found  by  fitting  a  smooth  curve  through 
the  experimental  values  obtained.  As  was  mentioned  above,  this  error 
analysis  was  performed  on  many  experimental  values  and,  from  observation 
of  these  results,  it  seems  reasonable  to  attribute  a  fractional  error  of 
cr  /  =  0.05  (or  5  percent)  to  all  unattenuated  gamma  dose  rates  re¬ 

sulting  from  the  sources  that  were  placed  on  the  roof  of  the  structure. 

F.  ERROR  IN  CONCRETE  THICKNESS 

•  The  standard  deviation  of  the  concrete  panels  used  in  this  series  of 
experiments  has  been  determined  to  be  about  1.3  percent^  of  the  nominal 
value  of  the  concrete  thickness.  To  determine  a  simple,  but  reasonably 
accurate  method  for  estimating  the  effect  of  the  concrete  thickness  error 
on  the  total  specific  dose  rate  within  the  structure,  certain  approximations 
were  made.  The  dose  rate  beyond  a  concrete  slab  of  thickness,  x  is  re¬ 
lated  approximately  to  the  unattenuated  dose  rate  by  the  equation: 


D(x)  =  D(0)e”^  (B-7) 

where: 

D(x)  =  attenuated  dose  rate  within  the  structure(R/hr/curie/ft2) 
D(0)  =  unattenuated  dose  rate  (R/hr/curie/ft2) 


Y.  Experimental  Calibration  and  Analysis  of  the  Radiation  Test 
Facility  at  the  Protective  Structures  Development  Center,  by 
C,  H.  McDonnell  et  al,  PSDC-TR-14  (September  1964) 
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p  =  linear  attenuation  coefficient'  ( (length) 
x  =  concrete  thickness  (units  of  length). 


Note  that  buildup  within  the  concrete  and  the  geometrical  effects  of 
the  structure  and  the  source  field  have  been  neglected.  This  is  a  fair  ass¬ 
umption  since  the  slope  of  the  attenuated  dose  rate  does  approximate  an 
exponential  of  this  type. 

The  total  differential  of  D(x)  is: 


AD<X>  -  •£$  AD  W  +  lir1  A* 


or 


AD(x)  =  e”^  X  AD(0)  -pD(0)  e"H  X  Ax 


This  is  equivalent  to 

AD(x)  _  AD(0) 

~m~  OW 


|1  Ax 


In  terms  of  fractional  error  this  becomes 


ffD(x) 

D W 


1/2 


(B— 8) 


where  px  equals  the  number  of  mean  free  paths  in  the  concrete  between  the 
source  and  the  detector.  The  quantity  cr^w  D(0)  is  the  some  quantity  as 
previously  evaluated  (see  Eq.  (B-5).  The  standard  deviation  of  concrete 
thickness  expressed  as  a  fractional  error  *r/  x  which,  as  mentioned  above, 
was  found  to  be  0,013. 

Thus  if  a  conservatively  high  value  of  aD(x)/  =  0.05  is  assumed, 
Eq.  (B-8)  becomes: 


gD(*) 

%) 


r, 

L 


0.0025  + 


(B-9) 


=  roof  reduction  factor 


Therefore 


specific  dose  rate  within  the  structure  (R/hr/burie/ft^) 
infinite  field  dose  rate  (R/hr/fcurie/ft^) 


=  W  iCM  +  Ic 


or  equivalently 


AD(x) 

~wr 


A  D0 

~TT 


In  terms  of  fraction  '  error: 


'*  -  (’OH  f  ,  (’Do 

R - |lW/  d0 


2  1/2 


crQ  /  Dq  is  equal  to  0.03685  and  crDx/  Dx  is  dependent  on  the  thickness 
or  concrete  between  source  and  detector.  Thus: 


R 

T~ 


_  faD(x)  \ 

"  ter  J 


+  0.001358 


The  total  fractional  error  inherent  in  lhe  roof  reduction  factor  as 
a  function  of  concrete  thickness  is  presented  below: 


on  and  Analysis  of  the  Radiation  Test  Facility 
at  the  Protective  Structures  Development  Center,  by  C.  H.  McDonnell, 
etal,  PSDC-TR-14  (September  1564). 

1,  An  Experimental  Evaluation  of  Roof  Reduction  Factors,  by  C.  H, 
McDonnell,  J.  Velletri,  PSDC-TR-16  (May  1566). 
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In  the  experiments  performed,  the  total  concrete  thicknesses  varied 
from  4  to  24  inches,  depending  on  which  structure  was  under  investigation 
and  on  which  floor  the  detector  in  question  was  located.  The  standard  de¬ 
viations,  expressed  as  fractional  errors  in  the  specific  dose  for  all  of  the 
concrete  thicknesses  used,  are  calculated  below : 


X 

(in.) 

px 

ir  °-ff 

*0  (x) 

WT~ 

0 

0 

0 

0.0$0  “ 

4 

1.34 

0. 000303 

0. 053 

8 

2.68 

0. 001214 

0.061 v 

12 

4. 02 

0.002731 

0.072 

16 

5.36 

0.0C4855 

0.086 

24 

8. 04 

_ 

0.010924 

0,116 

. 

Thus,  while  there  can  be  no  single  e.ror  attributed  to  the  specific 
dose  rate  in  the  structure  caused  by  roof  sources,  a  single  error  may  be  as¬ 
sociated  with  each  floor  thickness  under  consideration. 

G.  APPLICATION  OF  ERROR  ANALYSIS  TO  ROOF  REDUCTION  FACTORS 

Most  of  the  results  of  the  experimental  roof  programs  in  this  report- 
have  been  presented  as  roof  reduction  factors.  The  roof  reduction  factor  is 
defined  as  the  ratio  of  the  dose  rate  at  any  point  in  the  structure  caused  by 
contaminated  field  on  the  roof  to  the  dose  rate  at  a  point  3  feet  above  the 
ground  caused  by  an  infinite  plane  source  of  contamination.  The  infinite 
field  value  DQ  has  been  measured  in  previous  experiments^  to  have  a  value 
of  464  R/hr/curie/ft^.  There  is  a  fractional  error  of  crDo/D0  =  0. 03685 
associated  with  this  value.  The  evaluation  of  aDo/D0  can  be  found  in 
Reference  1,  with  a  slight  change  in  the  final  value  caused  by  a  later  cal¬ 
culation  of  the  cal  ibration  constant  error. 

The  calculation  of  aDo/D{)  includes  the  error  associated  with  a 
ground  roughness  correction  which  is  necessary  because  the  experimental 
measurement  was  of  necessity  conducted  on  other  than  a  mathematically 
smooth  plane.  The  total  fractional  error  in  the  roof  reduction  factor  may 
be  found  by  starting  with  the  definition  of  such  a  roof  reduction  factor: 

Rf=^i  (B-10) 
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X 

(in.) 

“dm 7  d(x> 

* 

VRf 

0 

0.U5 

0.062 

4 

0. 053 

0.065 

8 

0.061 

0.071 

12 

0. 072 

0.081 

16 

0.086 

0.094 

«$ 
_ J. 

0.116 

0.122 

The  standard  deviation  in  the  roof  reduction  factor  is  fc  jnd  by 
multiplying  the  reduction  factor  at  the  point  in  question  by  the  appropriate 

aRj/Rf 

H.  EXPERIMENTAL  VERIFICATION  OF  ERR  )R  ANALYSIS 

The  error  analysis  presented  above  is  based  on  estimates  of  errors 
or  uncertainties  in  the  instruments  used  and  other  parameters  associated 
with  determining  an  experimental  value  of  the  gamma  dose  rate.  The  ex¬ 
periments  conducted  at  the  Protective  Structures  Development  Center  are 
repeated  many  times  for  each  structural  configuration  so  that  a  variety  of 
points  are  available  to  yield  a  mean  value  of  the  dose  rate  at  any  particu¬ 
lar  position.  Figure  B-l  shows  typical  examples  of  the  various  experimental 
points  taken  on  two  floors  of  a  structure  together  with  the  mean  value  of 
the  dose  rate  and  the  limits  of  the  error  band  predicted  by  the  above  error 
analysis  in  attenuated  dose  rate.  The  derived  band  of  error  or  uncertainty 
seems  to  be  quite  good  in  predicting  experimental  error. 
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V.  ical  Position  (ft) 


FIG.  B-l  -  Typical  Repeated  Data  and  Error  Band 


APPENDIX  C 


EXPERIMENTAL  DATA 


Presentation  of  results  as  roof  reduction  factors  assumes  a  measured 
value  of  the  dose  encountered  at  a  point  3  feet  in  the  air  above  an  infinite, 
uniformly  contaminated  plane  of  radiation.  At  the  Protective  Structures 
Development  Center,  this  value  has  been  measured  to  be  464  R/hr/curie/ft^. 
A  detailed  dsscription  of  the  measurement  of  this  specific  dose  rate  is  pre¬ 
sented  in  Reference  1. 

Measurements  of  the  gamma  dos.  vvithin  the  structure  for  the  series 
of  experiments  analyzed  in  this  report  were  made  with  Victoreen  Model  362 
(200  mr),  Mode!  229  (10  mr),  and  Model  208  (1  mr),  nan-direct  reading 
ionization  chambers  (dosimeters)  together  with  o  Technical  Operation  Model 
556  charger-reader.  Selection  of  a  particular  dosimeter  for  any  specific 
experiment  was  based  on  the  dose  expected  at  the  detector  location  from  a 
particular  source,  source  field,  exposure  time,  and  concrete  thickness  be¬ 
tween  source  and  detector. 

The  dosimeters  used  and  the  charger-reader  were  calibrated  before 
use  in  the  experiments  against  a  gamma  source  of  known  strength  and  Vic- 
Hreen  R  meters  calibrated  by  the  National  Bureau  of  Standards.  The  only 
dosimeters  used  in  the  experiments  were  those  which  responded  to  within 
+  2  percent  of  the  known  dose.  At  various  times  during  the  experiments, 
a  secondary  calibration  was  performed  on  the  instruments  to  check  the  a- 
greement  of  the  dosimeters  with  each  other. 

Two  Co-60  gamma  sources  were  utilized  in  this  set  of  experiments. 
These  sources  had  values  of  4. 14  and  41. 4  curies  on  January  2,  1 966. 

The  dosimeter  readings  were  all  normalized  to  specific  dose  rates 
in  roentgens  per  hour  for  a  source  density  of  1  curie  of  Co-60  per  square 
foot  of  source  plane  area.  This  is  accomplished  by  means  of  the  following 
equation: 


D 


s 


MAT 

10^  p  C  5  t 


(C-l) 


where: 

Ds  =  specific  dose  rate  R/hr/curie/ft^ 
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/v\  =  Scale  reading  on  the  dosimeter  charger-reader  (pA) 

2 

A  =  area  of  simulated  source  field  (ft  ) 

T  =  temperature  (°R) 

p  =  barometric  pressure  (in.  Hg) 

pa 

c  =  charger-reader  calibration  constant 

t  =  exposure  time  (hours) 

s  =  source  strength  (curies) 

The  specific  dose  rates  measured  in  the  experiments  are  tabulated 
in  this  Appendix.  Table  C-l  shows  the  centerline  values  as  a  function  of 
source  detector  distance  for  both  the  full  roof  source  and  the  11.5  foot 
diameter  circle  for  the  four  and  eight  inch  floor  cases.  Tables  C-2  and  C-3 
present  similar  values  for  the  off-center  positions. 


°R 

in-Hg 


TABLE  C-l 


SPECIFIC  DOSE  RATES,  FULL  ROOF  SOURCE,  AND  DISK  SOURCE 

(R/hr/curIe/ft2) 

Center!  ho  Positions 


_ 48,6  psf  Roof  and  floors  97, 2  psf  Roof  and  floors _ 

Distance  [^stance 

from  source  from  source 

(ft)  Full  Roof  11. 5'd!a. circle  (ft)  Full  Roof  11.5'clule 

■  ■  ■■■>  ■  1  1  m ■  ■  "■■■■■■  ■  ■  ■  -  . . . .  ■  — —  —  ..I.  i  ■■  i  ■■  — — — . 


1.3 

16.4 

2 

20.9 

3 

21.4 

4 

21.2 

5 

21.3 

6 

20.1 

7 

19.0 

8 

18.4 

9 

17.5 

10 

16.7 

11 

16.0 

14 

3.18 

15 

3.30 

16 

3.25 

17 

3.13 

18 

2.97 

19 

2.77 

20 

2.66 

21 

2.49 

22 

2.38 

23 

2.27 

25.33 

0.364 

26.33 

0.442 

27.33 

0.464 

28.33 

0.467 

29.33 

0.453 

30.33 

0.435 

31.33 

0.413 

32.33 

0.402 

33.33 

0.390 

34.33 

0.408 

35.33 

0.359 

12.9 

16.7 

2 

15.8 

3 

14.5 

4 

11.3 

5 

9.57 

6 

7.90. 

7 

6.60 

8 

5.83 

9 

4.91 

10 

4.45 

11 

0.730 

14 

0.760 

15 

0.705 

16 

0.660 

17 

0. 580 

18 

0.550 

19 

0.509 

20 

0.480 

21 

0.447 

22 

0.411 

23 

0.0712 

0.0780 

26.67 

0.0770 

27.67 

0.0800 

28.67 

0.0722 

29.67 

0.0686 

30.67 

0,0645 

31.67 

0.0619 

32.67 

0.0608 

33.67 

0.0568 

0,0542 

34.67 

4,50 

4,2 

5.20 

4.25 

5.30 

4.04 

5.30 

3.52 

5.23 

2.96 

5.03 

2.54 

4.92 

2,20 

4.76 

1.90 

4.55 

1.66 

4.54 

1.46 

0.206 

0.0812 

0.236 

0.0835 

0.246 

0.0830 

0.246 

0.0760 

0. 238 

0.0675 

0.227 

0.0625 

0.215 

0.0560 

0.205 

0.0516 

0.197 

0.0487 

0.192 

0.0440 

0.0112 

0.0118 

0.0118 

0.0118 

0.0109 

0.0106 

0.0103 

0.0100 

0.00916 
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TABLE  e-a 


SPECIFIC  DOSE  RATES,  FULL  ROOF  SOURCE 
.  /hr/curie/ft^) 


Off  Center  Positions 
48.6-psf  floors 


Position  *  Distance  from  source  (ft) 


X,Y  ~ 

w 

i\ 

18 

15 

- ^ - 

+10, +15 

0.220 

0.244 

0.298 

1.26 

1.41 

1.70 

9.18 

U.7 

18.1 

-10, +15 

0.234 

0.256 

0.3T8 

1.29 

1.47 

1.71 

9.12 

11.2 

16.3 

-10,-15 

0.225 

0.255 

0.309 

1.32 

U41 

1.70 

8.76 

10.7 

16.0 

+10,-15 

0.236 

0.252 

0.318 

1.52 

1.50 

1.74 

9.47 

11.2 

16.8 

46, +9 

0.329 

0.375 

0.465 

1.99 

2.33 

2.81 

15.2 

17.7 

20.7 

-6, +9 

0.256 

0.365 

0.440 

2.01 

2.37 

2.76 

14.8 

17.9 

21.3 

-6,-9 

0.320 

0.364 

0.452 

1.99 

2.37 

2.87 

14.9 

17.9 

21.2 

•»6,  -9 

0.326 

0.375 

0.460 

2.02 

2.37 

2.88 

15.0 

17.6 

21.3 

0,+15 

0.310 

0.341 

0.353 

1.71 

2.06 

2.31 

12.9 

15.7 

20.2 

0,-15 

0.302 

0.333 

0.352 

1.68 

2.02 

2.35 

12.5 

15.3 

19.5 

0,+9 

0.348 

0.395 

2.16 

2.46 

2.79 

16.5 

18.9 

21.4 

0,-9 

0.348 

0.393 

— — 

— — 

16.1 

— 

21.0 

46,0 

0.372 

0.432 

0.515 

2.23 

2,72 

3.13 

15.9 

18.4 

21.3 

-6,0 

0.366 

0.423 

0.504 

2,19 

2.76 

3.32 

16.1 

18.7 

21.8 

+10,0 

0.311 

0.348 

0.425 

1.77 

2.11 

2.47 

12.5 

14.2 

18.6 

-10,0 

0.318 

0.360 

0.447 

1.83 

2.10 

2.55 

12.5 

14.3 

18.9 

♦See  Figure  3.4 
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tabu  C-3 

SPECIFIC  DOSE  RATES,  FULL  ROOF  SOURCE 

(R/hr/curie/ft2) 

Off  Center  Positions 
97. 2-psf  floott 


Position  • 

Distance  from  Source  (ft) 

X,  Y 

21 

T8 

15 

9 

6 

3 

+10, +15 

_ 

0.108 

0.143 

2.65 

3o40 

5.20 

-10, +15 

0.090 

0.102 

0,131 

2.4* 

341 

4.63 

-10,-15 

0.102 

0.102 

0.137 

2.40 

3.11 

4.50 

+10,-15 

0.0% 

0.102 

0.137 

2.38 

3.04 

4.45 

46, +9 

0.167 

0.195 

0.260 

4.00 

4.60 

5.40 

-6, +9 

0. 165 

0.193 

0.231 

3.% 

4.50 

5.25 

-6,  -9 

0.163 

0.193 

0.236 

3.97 

4.61 

5.40 

46,-9 

0.159 

0.19C 

0.225 

4.04 

4.63 

5.60 

0,+15 

0.154 

0.166 

0.190 

3.58 

4.10 

5.20 

0,-15 

0.143 

0.161 

0.178 

3.40 

4.02 

4.95 

0,+9 

0.178 

0.207 

0.231 

4.24 

4.80 

5.15 

0,-9 

0.178 

0.207 

0.231 

4.24 

4,62 

5.00 

46,0 

0.184 

0.212 

0.242 

4.12 

4.77 

5.41 

-6,0 

0.178 

0.212 

0.254 

4.20 

4.90 

5.50 

+10,0 

0.143 

0.161 

0.178 

3.42 

4.00 

4.% 

-10,0 

0.143 

0.172 

0.1% 

3.34 

3.82 

5.03 

•See  Figure  3.4 
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